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Abstract 

Ion Mobility Spectrometry (IMS) is a well-known, sensitive and rapid technique 
to detect dangerous organic compounds. We propose a system in which a crown 
type discharge generates a ionic flux that is swept towards an array of collectors 
by a transverse electric field. The ions are separated as they enter the cell ac- 
cording to their mobility. Thus, the distribution of the charge collected at the 
detector assembly constitutes a fingerprint for each organic compound. Simula- 
tions of our cell and experiments were performed for small amount of Acetone, 
Ethanol and Toluene. The dependence on the cell parameters of the current 
and charge versus time of flight was analyzed. Our simulation reproduces only 
qualitatively the experimental results. However, a PGA statistical analysis of 
the simulated results shows that such a fingerprint can be used to identify those 
compounds with certainty. 



1. Introduction 

The technique for identification of volatile organic compounds by differences 
of ion mobility is one of the most efficient techniques applied for detection of 
explosives, drugs and industrial toxic compounds The detection limit of this 
technique is found to be of the order of 10~® (ppb). It has a very rapid response 
time, which makes it very attractive for situations of high control demands. The 
ions of the molecule to be identified acquire a velocity under the influence of an 
external electric field E. At the lowest order 
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VE = KE, (1) 

where K is the mobiUty. If the field is intense enough it has been reported 0] 
that one should go two orders higher (cubic) in the field dependence. A way 
to understand this result is by considering a power series expansion on E. As- 
suming isotropy of K and applying symmetries on spatial transformations, the 
vectorial properties of ve and E determines that only odd terms in such expan- 
sion are possibles. 

Besides the interaction of the ions with the field E, they collides with 
molecules of the carrier gas that fiow through the detection system under a 
laminar flux. The explicit dependence of K on the pressure due to the carrier 
gas, the temperature, the ion charge, its mass and the cross section can be 
obtained from general considerations when the field E is not very intense 
Working on the center of mass framework to describe the collisions between 
two bodies, and the reduced mass rrir — Mfn/{M -\- m) where m is the mass 
of ions and M the mass of the molecules of the gas, the variation of the linear 
momentum of ion-gas system is due to the external impulse of force, i.e. qEr, 
being q the ion charge and t the characteristic time between collisions. It can 
thus be established that 



Ve — qEr/mr. (2) 

Furthermore r can be related to a collision cross section of the ion-gas system fl 
and the density of this system, which is approximately equal to the gas density 
rj. Considering the characteristic length I = vt and the average rate v of the 
reduced system ion-gas between collisions, it is estimated that in the volume l^l 
a collision should typically occur. Hence -q « l/{lVt) and Eq. ([2]) can be written 

qE 1 

VE ^ pp. (3) 

One can approximate v « {iksT /mrY/"^ because the kinetic energy of the ion- 
gas system is mainly due to thermal agitation of value /cbT/2 for each degree 
of freedom. Then 

qE ( m + M \ ^''^ 
~ rjQ \mM3kBT ) ' ^ ' 

being ks the Boltzman constant and T the absolute temperature. In Eq. ([4]) 
there is a factor 3~^/^ that should be replaced by 3(27r)^/^/16 when more ad- 
vanced molecular calculations are considered 0]. 

We include such derivation here because it contains the right dependence 
with the physical parameters of interest and its derivation is straightforward. 
In the calculations presented in this work we employ instead the mobility usually 
reported in the literature 0, Q , which is obtained from 

3q ( 2n{m + M) \''^ 
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This mobility is standardized at normal Pq pressures and Tq temperatures 
through 

Ko^KiTo/T,P/Po), (6) 

To — 273 °K and Pq depends on the system of unities used to refer to 1013 HPa. 
Kq is known as the reduced mobility. 

A decade ago the identification of species by separation techniques l6j] was 
based on the differences in the reduced mobilities given by Eq.®. Since then 
substantial modifications have been made in the design and size of sensor de- 
vices. Using asymmetric electric fields in the range of radio-frequency and by 
micro-machining the detection cell, a field asymmetric or FA-IMS variant was 
proposed Identifications of up to 0.1 ppm of Toluene and also a resolution 
threshold of 60 ppb for a mix of gases of Acetone-Benzene and Acetone- Toluene 
were reported. It has been recently proposed 0] a combination of mass spec- 
troscopy and ion mobility to improve the transmittance in an Ion Mass Mobility 
Spectroscopy through a focus cell by means of periodical electrostatic fields dis- 
tributed on a very long arrangement of electrode rings (Im of length). It has 
also been proposed recently an IMS type cell by swept-field scanning ^ in which 
a plane and parallel detectors system produces an electric field perpendicular 
to the flux of the ions being detected. The operating principle of this device for 
mobility identification is based on the differences of the average path described 
by the ions within the detection cell. 

F = 



Vd 

\V = Vo 



Figure 1; Average paths within the detection cell. Di, Dj represent two of the seven detectors 
contained in the cell. The detectors are earthed and the anode is kept at Vq potential. The 
left path corresponds to an ion with higher mobility with respect to the right path. 

In Fig. [T] an example of two ion paths corresponding to ions with different 
mobility is given. Due to the velocity of the carrier fiux, the ions reach a de- 
tection zone with velocity Vd parallel to the detectors plane. By the action of 
the transversal field E, the ions acquire, according to Eq. ([1]), a vertical compo- 
nent velocity ve proportional to the mobility K. By adjusting the geometrical 
parameters of the cell, the flux of the carrier gas and the potentials among elec- 
trodes, it is possible to obtain a typical ion distribution for each compound on 
the Di and Dj detectors, where i ^ j = 1 . . .n are labels. This is characterized 
by a greater accumulation of charges on the detectors, whose positions measured 
from the gas input in the cell are more distant for ions with less mobility, as is 
suggested schematically in Fig. [TJ 

A modification in the swept-field scanning detection was recently proposed. 
It includes an ion-focusing just before the introduction of ions in the detection 
cell [13] ■ In this variant, the flux of both the carrier gas and the ionized gas are 




3 



separated through geometrical constrictions. This increase the precision of the 
positions reached by the ions in the detection cell. In that sense, the aim of such 
a proposal is not meant to obtain a fingerprint for each compound. Indeed it 
proposes that the identification shall be determined by the average time taken 
by each species to reach a specific detector [l^ . Although the ion time of flight 
has been analyzed for different values of carrier gas flux and potentials between 
electrodes in order to evaluate the focalization conditions, it has been found in 
this work that identification of species by statistical analysis of their traces on 
detectors is more accurate if few restrictions are imposed by operational condi- 
tions, i.e., work pressure and non-radioactive ionization method. Furthermore, 
due to the fact of a very high difference of electric potential generated by the 
crown effect, fragments of the molecule to identify are produced, each one with 
a different mobility. In our work we have found that this effect is consistent 
with the formation of a fingerprint for each species at the moment of measuring 
the total charge deposited on the detector assembly. 

Basically, the main idea behind a fingerprint identification is the following: 
A system is designed in such a way that traces of given chemical components are 
well defined under specific operating conditions. With this purpose, numerical 
simulations are first performed, which allow us to identify these specific condi- 
tions. The algorithm used to simulate the ion detection system has proven to 
be sufficiently robust for the design of ion mobility detection cells Fur- 
thermore in our work, the design of an experimental prototype is analyzed for 
obtaining traces of Ethanol, Acetone and Toluene. The statistical study of 
these traces was performed by the principal components analysis [13] , known as 
PCA (acronyms from Principal Components Analysis), which was used earlier 
on similar identifications within the design of an electronic sniffer [l^ . 

We present our work in the following way. In Section [2] we describe the ex- 
perimental configuration of the prototype device and the procedure we followed 
for the measurement of signals in the identification of three different chemical 
species. In Section [3] we present a brief description of the algorithm used for 
the simulation of ion paths. We divided Section [5] in four parts in order to 
show the experimental results in Subsection 14.11 time of flight simulations in 
Subsection 14.21 and collected charges in Subsection 14.31 In Subsection 14.41 we 
present the analysis of principal components and Section [5] is devoted to the 
main conclusions. 

2. Experimental device 

An experimental device has been designed as indicated in Fig. [2l It consists 
of three principal parts: i) the carrier gas flux inlet, ii) the ionization chamber 
and iii) the detection cell. The sample is introduced as a gas in the ionization 
chamber together with the carrier gas where, through a type of crown discharge, 
is ionized. The pointer shaped electrode, shown in the central part of Fig. [21 
is a stainless steel needle of 0,8 mm of diameter fixed inside a teflon cylin- 
der. An earthed copper tube of 8 mm of diameter closes the discharge circuit 
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Figure 2: Lateral section of the experimental device. The upper part corresponds to the 
carrier gas inlet. In the central part ionization is produced, and the detection cell is on the 
right. 

of the pulse. The ionization discharge is recorded with an oscilloscope (Tek- 
tronix P6015). The needle is concentric to the copper tube; hence the distance 
between electrodes is approximately 4 mm. A stainless steel mesh is placed 
between the ionization chamber and the detection cell in order to screen the 
high field produced during the ionization discharge within that last zone. The 
detection arrangement is similar to that proposed by Soils et al. Q ; the field E 
is transversal to the ion flux determined by the carrier gas. Our ion detection 
device is an assembly of 7 earthed detectors placed in a plane as shown in Fig. [2] 
and connected to a constant potential which is fixed in the range between 10 and 
lOOV. The dimensions of the detection cell are 17mm x 11mm x 2mm. A dry 
carrier gas of N2 is used in order to avoid humidity problems in the air (Indura 
99,998%, H2O < 3mg/ml). The gas flux is controlled by a flux measuring device: 
Brooks model 6-658. The testing organic substances are absorbed in a piece of 
cotton placed in a glass tramp. A flux of N2 gas passes through the tramp drag- 
ging the steam of the sample at room temperature and directed to the ionization 
zone. The ions generated are dragged by the same flux of N2 gas towards the 
detection cell. Each electrode is connected to a resistance of 3.2 Mfi and the 
electrical signal generated by the ions is measured in a oscilloscope Instck GDS 
810S synchronized with the pulse of the discharge of the ionization needle. All 
the system is controlled with a PC and due to the high signal/noise relationship, 
only 8 signals are necessary to average the response of each detector. With this 
system it is possible to study: Absolute Ethanol (Merck, pro-analysis, 99.8%), 
Toluene (Mallinckrodt, analytical, 99.8%) and Acetone (Sintorgan, analytical, 
99.5%). All of them were used without any previous treatment. 
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3. Numerical simulations 

The ion paths in the detection ceU were simulated with the SIMION 8.0 
code. The SDS module was implemented. This module combines the dynamics 
of a particle in a viscous medium and its interaction with the carrier gas at room 
pressure We give here a brief description of such a simulation. In order to 
define an ion path, the equation of motion is considered, in which acceleration 
has a viscosity component due to velocity in Eq.(IT]), and another Coulombic 
component due to the electrostatic potential of the cell. Furthermore when the 
SDS module is implemented, a diffusion model is incorporated. The statistics 
of the distance covered after 10^ collisions in 10^ ion samples was applied. From 
the ion-gas mass relationship data of 1, 10, 100, 1000 and 10000, the distance 
corresponding to intermediate values of mass were calculated by interpolation. 
Such distance is used to modelling the size of the jump in a step of the path. The 
direction within the 3D space of the cell should be calculated with a uniform 
distribution in order to modelling the expected local isotropy. In this sense, an 
appropriate sample space should be chosen. For instance, if a point within a cube 
is randomly taken in order to obtain a random direction, the appropriate sample 
space is defined only by those points that belong to the sphere inscribed in such 
a cube [iJl- Once the path is defined, the algorithm reports the positions of the 
ion's impact. We have instrumented a couple of codes written in interpretative 
language (Perl) to process the data accompanying each ion path. With these 
codes we obtain the charge collected as a function of time. Through a numerical 
derivation of the collected charge, the current in each detector as a function of 
the time of fiight was obtained. 

The SIMION 8.0/SDS algorithm ^] allows us to estimate (from the mass 
data) or to take as a predeterminated information, the mobilities of the ions 
whose paths will be modeled. Two types of ion path simulations were consid- 
ered: one without fragmentation and the other with fragmentation. This means 
that, for each species studied, ions of the whole molecules were considered first, 
and secondly, the fragments of such molecules according to their relative abun- 
dances [15] being the least of them higher than 10%. Fig. [3] shows the surface 
potential energy view obtained with Simion 8.0 for the simulation of paths con- 
sidering only primary ions (without fragmentation) of Toluene, Acetone and 
Ethanol. Results indicate that the ions of Acetone and Ethanol with greater 
mobility are collected preferentially in the detectors closer to the capture cell 
inlet (see also Fig. [2]). 

4. Results 

In this section, experimental results and numerical simulations of the device 
proposed here for the detection cell by swept-field scanning are presented. Flight 
times are analyzed as well as the collected charges in detectors for different 
combination of carrier gas flux and detector potentials. The analysis of the 
principal components for one of the selected combinations is also presented. 



6 



Figure 3: Surface potential energy view of the ion path simulation of a Toluene (black) - 
Acetone (light gray) - Ethanol (gray) mixture obtained by Simion 8.0. The Ions come from 
the rigth top side of figure. The top part of the traces correspond to the potential of needle 
discharge (5000V). The carrier flux with velocity 30 m/s allows the flow of the ions go into 
detector zone where potential increments are indicated by the surface potential level. In the 
left bottom side of figure ions impacts to detector at OV at the lowest part of the surface 
potential. The anode is opposite to it in the upper part of the surface potential at 50V. 



The experimental results are given in the first subsection, and in the following 
subsections the results of simulations and numerical calculations are shown. 

4.I. Experimental 

Under normal conditions of pressure and temperature different tests were 
performed as described in Section [2] for Acetone, Ethanol and Toluene at satu- 
rated steam pressure using the device shown in Fig. [51 The location of the detec- 
tors in the detection cell is shown in Fig. [4l Ion gas moves from left to right as 



















Figure 4: Diagram that shows the numbering for the array of detectores. Ions flux goes into 
detection cell from the left as indicated by the arrow. 

indicated by the arrow. The charge collected during the detection process was 
obtained for each detector. In Fig.[5]these data for each substance are displayed. 
According to Eq.(I2), the differences in mobility allows to correlate time with 
position taken by an ion which arrives on a detector. In order to explain exper- 
imental results, we need to consider that the mobility of Toluene {Kq = 1, 87) 
is lower than that of Ethanol {Kq = 2, 06) and Acetone (i^o = 2, 13) (j, |j 
and that with the configuration proposed, the electric field is perpendicular to 
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Figure 5: Collected charge measured in each detector when Acetone (light gray), Ethanol 
(gray) and Toluene (dark gray) samples are tested separately. 



the direction of the carrier gas (see Fig. [T]). Hence, Toluene ions should travel 
longer distances without being deviated compared with Ethanol and Acetone 
ions. For these last ones, the scarce relative difference of 3.5% reported be- 
tween the Acetone (larger) and Ethanol (smaller) mobilities 0, rnight yields 
to similar path deflection and nearly the same distance of impact on detectors 
produced by the swept-field scanning. In general, it is expected that lighter ions 
have larger mobilities as they have smaller inertial mass. Although this rule ap- 
plies for a large number of molecules Q there are some exceptions where the 
opposite is true. Such exceptions can be explained in terms of a fragmentation 
process of the primary ion to lighter secondary ions and, therefore, with greater 
mobility. Due to the type of ionization proposed in this work, we suppose that 
the fragmentation of Acetone and Ethanol produce secondary ions with highly 
relative abundance [l6j . Consequently, each species aquire a different mobility 
compared to that expected in terms of primary ions mass alone. 



4-2. Times of Flight 

The velocity of the carrier gas and the difference of electric potential within 
the cell can be conveniently combined in order to modulate the amount of ions 
that goes to a given place of the experimental device. We shall analyse here 
the time of flight for each ion and the average rate of the collected charge 
in each detector as a function of both the velocity of the carrier gas and the 
difference of electric potential. In Fig. [6] we display the normalized current 
generated in each detector by the ions of Toluene when the velocity of the 
carrier gas is 3, 9 and 30 m/s and the differences of electric potential in each 
detector cell are 15 and 150 V. A pattern that spreads the current over time 
of flights when both the velocities Vg and potentials diminishes is observed. 
This can be explained considering that the contribution to the ion's energy 
from the interaction between the ions and the carrier gas do compete with 
the electrical counterpart. Increasing the electric potential differences of the 
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Figure 6: Normalized I/Imax current intensity vs. time of flight (TOF) for each detector Di 
with i = 1 ... 7 ordered by pairs as shown in the graph labels. For Vg = 3 and 9 m/s, the 
difference of potential is A = 15 V. 



cell will give a larger kinetic energy contributtion to the ions compared with the 
contribution from the Brownian motion. Thus, such electric potential increment 
clearly produce a lesser diffuse and so more defined and streched ion paths. For 
that reason, in Fig. [6] it is observed that, as both the velocity of the carrier gas 
and the potential decrease, the average rate of ions that arrive at the same time 
also decrease. This means that in such case the ions arrival are expanded over 
the whole detection zone. 

4-3. Collected charges in detectors 

Counting the amount of collected charge in each detector offers an alterna- 
tive way to get the times of flight. Their results are also in agreement with 
experimental data. We present here the measurement of the collected charge 
when the velocity of the carrier gas is 30 m/s and the difference of potential 
in the cell is 50 V. ft is verified that the times of flight are larger for ions that 
arrive on detectors which are farther to cell inlet. For instance, in Fig. [7] the 
same graph for detectors Dg and Dj is shown for the collected charge Q di- 
vided by total charge Qt of ions when Toluene, Acetone and Ethanol are tested 
separately. 

When mobility is estimated from ion mass, as should be if we consider the 
case of fragmentation process during the ions generation, the charge collected 
in each detector as a function of time of flight follows the simple rule expected 
by inertial mass considerations. This can be observed in Fig. [T] The collected 
charge for Ethanol is greater than that of Acetone for times of flight between 280 
and 350fis corresponding to the arrival of ions at the detector Dq. However, if we 
consider the mobilities reported in the literature we could expect that this result 
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Figure 7: Normalized total charge Q/Qt for detectors Dg (left) and Dy (right) for each 
species as a function of time of flight when primary ions were considered for Toluene (full 
lines), Ethanol (dotted lines) and Acetone (dashed lines). 



should be the opposite one. For that reason, we performed our simulatfon also 
considering groups of secondary ions. We supposed that when Acetone, Ethanol 
and Toluene were ionized they would form groups of ions from the fragmentation 
of the primary ion. We considered the ions with a relative abundance greater 
than 10%. Furthermore, using secondary ions, we found out a large defocusing 
and a low temporal definition due to a wide charge distribution on the detectors 
assembly. Then as an improvement we considered more appropriate to represent 
the collected charges in each detector taking care also of the secondary ions for 
each species. In Fig.|8]we display such results. In this case a charge distribution 
closer to that expected in terms of Acetone ion mobilities is obtained. However, 
a larger amount of Ethanol ions than that of Acetone ions in detectors {Di,Dr,) 
is still observed. Such detectors are the first pair from the entrance to the cell 
of the carrier gas. The second pair of detectors is {D2,Dq). For this second pair, 
the method of secondary ions seems to be more consistent with the expected 
result in terms of known mobilities. The detector pair that follows is (D^jDr). 
In this last case, the result that arise when considering fragmentation is again 
opposite to that which would be expected from mobilities. 

4-. 4- Principal components analysis 

Based on the previous background on simulation and experimental results 
we are be able do another step forward. Our working assumption is now that 
once the cell geometry and the operational conditions of our device are given 
(the carrier gas flux and the cell potential values) the charge collected in each 
detector has its own statistics for each unknown species. Specifically, the analy- 
sis of the principal components applied to our device allows us to determine the 
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Figure 8: Normalized collected charge Q/Qt in each detector obtained from the fragmentation 
product of the primary ion with relative abundance greater than 10% for Acetone (light gray), 
Ethanol (gray) and Toluene (dark gray). Over 6000 ion paths are averaged for each species 
considering carrier gas velocities between 20 and 30 m/s. Numeration and configuration of 
detectors is the same as that of the experimental device. 



fingerprint of eacli species corresponding to our experimental configuration. We 
associate a random variable X* (i = 1 ... 7) to the collected charge Q in each 
detector, normalized with the total charge Qt and distributed on all the detec- 
tors. Thus, X'^ are real numbers in the interval [0, 1]. Each measurement of the 
collected charge in a detector corresponds to an event of the random variable X' . 
Therefore, if we perform TV measurements for each variable, we can construct 
the vector X* = . . . , whose coordinates are the result of these 

measurements. The notation (. . .)^ means the transpose of (. . .). From the 
definition of the mean value = ^/Nj^i ^-nd experimental covariance [l3|, 
a real symmetric and covariant matrix C is obtained, whose dimension is equal 
to the number of detectors, i.e. 7x7. Their elements are defined by 

EE Cov{X\X^) = —1- Y,iK. - X^){Xl - X^), (7) 

where m = 1 . . . N . We now introduce the similarity transformation 

D = P^CP, (8) 

where D is a diagonal matrix and P is an orthogonal matrix containing the 
eigenvectors of C ordered by columns. From the last matrix it is possible to 
define the transformation 

Y = P^X, (9) 
That enable us to change from random variables of detectors 

X= (10) 
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to a new set of random variables 



Y = (f^f^...,fT, (11) 

From this output a principal subgroup can be obtained 

Z = Q^X, (12) 

being Q the matrix formed by the first n < 7 eigenvectors of C ordered in 
decreasing order of its eigenvalues. The final results 

Z = (Z\...,Z"f, (13) 

contain the expected principal components. Different criteria were proposed for 
getting n [l2| , although all of them consider the minor subspace in which the 
greatest quantity of data is concentrated 

The results of Eq. pT|) can be understood in terms of a coordinate transfor- 
mation. Within the new coordinate system, each component is independent. 
This allows us to define an orthogonal space that represent the data obtained 
from measurements of collected charges for each species. The principal eigen- 
vectors used show in turn the specific combination of the original detectors that 
intervene in the output of the measurement. 
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Figure 9: The first two principal components for 5 measurements of tfie charge collected in 
each detector when using all of them. Ethanol (plus symbol), Acetone (cross symbol) and 
Toluene (asterisk symbol) averaged on 6000 runs each species. 



Fig. [5] shows the coordinates of the two first principal components for a set 
of 5 measurements of the collected charge in each detector for Ethanol, Acetone 
and Toluene ions and fragments of them. Each one of these measurements is 
similar to that represented in Fig. [8] and corresponds to the average values on 
a sample of 6000 ion paths, considering the carrier gas velocities between 20 
and 30 m/s. Only 2% of data remains outside the analyzed set by the two first 
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principal components, distributing 91% in the first and 7% in the second of 
them. Fig. [10] shows the coordinates in the two first principal components of 
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Figure 10: The first 2 principal components for 5 measurements of collected charge in each 
detector, but without detector number 4. Ethanol (plus symbol). Acetone (cross symbol) and 
Toluene (asterisk symbol) averaged on 6000 runs each species. 

measurements on the same system considered in Fig. ^ but with one variable 
less due to detector number 4 was not considered. In this case it is observed 
that 99% of data is represented by the first two components and 98% only by 
the first one. This is explained by the fact that being very scarce the collected 
charges on detector 4, as can be seen in Fig. |51 the diminution of the number 
of variables decreases the dimension of the total space considered with a very 
little alteration of the sample size. Therefore, the concentration of data in the 
first two principal components increase. The relevance of the PCA result is 
considerable regarding the separation of the three set of data for all species 
we have studied. The identification of an unknown sample is positive if its 
representation fits within any of the clearly non-overlapping region of our three 
set of data (see Figs. [9l and fTO|) . 

5. Conclusions 

Our simulation clearly shows that for the trace of ions, the carrier gas velocity 
and the cell potential are directly correlated. The analysis of fragmented ions 
partially explain the traces of Acetone and Ethanol when their mobilities are 
estimated from fragments mass. However, the integrated signals which were 
obtained experimentally show the same trend that our simulation. As shown in 
Fig. [5l when comparing the signals of detectors {Di,D^) with the subsequent 
pair {D2,Dq), it is observed that they maintain almost the same proportion 
between all substances for each detector, but the whole amount is increased 
in the second pair. In the third pair, (D^^Dj), the proportion of the trend 
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between substances almost remains the same, but decreasing the total amount 
of the set. Only in the most distant detector the result is different and shows the 
proportions between species following the expected results in terms of tabulated 
mobilities 

Since the simulations of the times of flight for each substance shows a char- 
acteristic distribution for both, charges and currents obtained through the de- 
tectors assembly, i.e. a fingerprint, the statistical correlation in the distribution 
of collected charges in the detectors using the Principal Component Analysis 
can be applied. 

For the experimental prototype proposed here a set of specific parameters 
(velocity of carrier gas and electric potential of the cell) to define Acetone, 
Ethanol and Toluene traces were optimized by numeric simulation. Afterwards 
a PCA statistical analysis of the collected ion charges in the set of detectors is 
used to define a non-overlapping region in a appropriated diagram for each of the 
above mentioned compounds. Thus, our procedure, which involves a training to 
develop the proper fringerprint of each compound, can be used in the designing 
of devices for the identification with certainty of unknown compounds, 
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